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ABSTRACT 

We present the discovery of two T dwarf binaries, 2MASS 1225—2739AB and 2MASS 1534—2952AB, 
identified in a sample of ten T dwarfs imaged with the Hubble Space Telescope Wide Field Planetary 
Camera 2. Companionship is established by the uniquely red F814W—F1042M colors of the binary 
components, caused by heavily pressure-broadened K I absorption centered at 7665 & 7699 A. The 
separations of the two binary systems are 0'/282±0'/005 and 0'.'065±0'/007, implying projected separations 
of 3.17±0.14 and 1.0±0.3 AU, respectively. These close separations are similar to those found in previous 
brown dwarf binary searches, and permit orbital mapping over the coming decade. 2MASS 1225—2739AB 
has a substantially fainter secondary, with AMi?8i4W = 1.59±0.04 and AMirio42M = 1.05±0.03; this 
system is likely composed of a T6 primary and T8 secondary with mass ratio 0.7-0.8. The observed binary 
fraction of our HST sample, is consistent with results obtained for late-M and L field dwarfs, and 

implies a bias-corrected binary fraction of for a > 1 AU and q > 0.4, significantly lower than the 

binary fractions of F-G and early-type M dwarf stars. Neither of the T binaries have separations a > 10 
AU, consistent with results from other brown dwarf binary searches. Using the statistical models of 
Weinberg, Shapiro, & Wasserman, we conclude that tidal disruption by passing stars or Giant Molecular 
Clouds, which limits the extent of wide stellar binaries, plays no role in eliminating wide brown dwarf 
binaries, implying either disruption very early in the formation process (ages <1 — 10 Myr) or a formation 
mechanism which precludes such systems. We find that the maximum binary separation in the brown 
dwarf regime appears to scale as a possible clue to the physical mechanism which restricts wide 

substellar systems. 

Subject headings: stars: binaries: visual — stars: formation — stars: fundamental parameters — stars: 

individual (2MASS J12171110-0311131, 2MASS J12255432-2739466, 2MASS 
J15344984—2952274) — stars: low mass, brown dwarfs 


1. INTRODUCTION 

T dwarfs are low-temperature {Teff ^ 1300-1500 K) 
brown dwarfs exhibiting distinct absorption bands of CH4 
in the near-infrared H- and K-bands (Kirkpatrick et al. 
1999; Burgasser et al. 2002d; Geballe et al. 2002). They 
are distinguished from warmer L-type (Kirkpatrick et al. 
1999; Martin et al. 1999) and M-type brown dwarfs by 
the presence of the CH4 bands, in addition to significant 
H2O and collision-induced (CIA) H2 absorption in the 
near-infrared. These molecular features, combined with 
pressure-broadened K I and Na I absorption at red-optical 
wavelengths (Tsuji, Ohnaka, & Aoki 1999; Burrows, Mar- 
ley, & Sharp 2000; Liebert et al. 2000), force the emergent 
spectral energy distributions of T dwarfs to peak around 
1 ^m. Since the discovery of the prototype of this class, 
Gliese 229B (Nakajima et al. 1995; Oppenheimer et al. 
1995), over 30 T dwarfs have been identified in wide-held 
surveys such as the Two Micron All Sky Survey (Burgasser 


et al. 1999, 2000c, 2003, 2002d, hereafter 2MASS; Skrut- 
skie et al. 1997) and the Sloan Digital Sky Survey (Strauss 
et al. 1999; Tsvetanov et al. 2000; Leggett et al. 2000; 
Geballe et al. 2002, hereafter SDSS; York et al. 2000), deep 
narrow-held surveys (Cuby et al. 1999; Liu et al. 2002; Za- 
patero Osorio et al. 2002), and as companions to nearby 
stars (Nakajima et al. 1995; Burgasser et al. 2000b; Els et 
al. 2001). 

Eor both stars and brown dwarfs, multiplicity is one 
of the fundamental properties that probes formation pro¬ 
cesses. Multiple systems are common among main se¬ 
quence stars, with roughly 60% of solar-type stellar sys¬ 
tems found to be comprised of two or more components 
(Abt & Levy 1976; Abt 1987; Duquennoy & Mayor 1991). 
In contrast, only 32-42% of held M dwarf systems are mul¬ 
tiple (Henry & McCarthy 1990; Fischer & Marcy 1992; 
Reid & Gizis 1997). Recent investigations of late-M and 
L held dwarf samples yield even smaller multiplicity frac¬ 
tions, only 20-30% (Koerner et al. 1999; Reid et al. 2001; 
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Close et al. 2002). Furthermore, while the separation dis¬ 
tribution of F-M stellar pairs appears to be broad (rang¬ 
ing from approximately 0.1 AU to 0.1 pc) and unimodal 
(peaking around 3-30 AU), all late-M and L dwarf binaries 
identified to date have apparent separations a < Oi'6 and 
projected separations Qproj < 15 AU (Martin, Brandner & 
Basri 1999; Koerner et al. 1999; Reid et al. 2001; Leggett 
et al. 2001; Close et al. 2002). Ejection models for brown 
dwarf formation (Reipurth & Clarke 2001; Bate, Bonnell, 
& Bromm 2002) have been proposed to explain this dearth 
of widely-separated, low-mass dwarf pairs, which may be 
the result of nature (i.e., inherent in the formation process 
itself) or nurture (i.e., due to dynamical scattering). 

In this article, we present imaging results for the first 
T dwarf binary search sample, undertaken using the Hub¬ 
ble Space Telescope (HST) Wide Field Planetary Camera 
2 (Biretta et al. 2002, hereafter WPFC2). In §2 we de¬ 
scribe the observations and image reduction techniques, 
and identify closely-separated sources near two of our T 
dwarf targets. In §3 we describe our photometric analysis, 
presenting aperture photometry and colors for our com¬ 
plete sample. The colors of the individual components of 
the 2MASS 1225-2739® and 2MASS 1534-2952 pairs con¬ 
firm their companionship. We discuss our PSF fitting in 
§4, by which we derive rigorous flux ratios and separations 
for the two T binaries, and quantify our search limits. In §5 
we discuss individual targets in detail, including a possible 
faint companion to 2MASS 1217—0311. In §6 we analyze 
binary statistics for field L and T dwarfs, and compare to 
stellar samples. Finally, we discuss our results in light of 
brown dwarf binary formation and destruction processes 
in §7, and conclude that the small separations of these sys¬ 
tems are not due to disruptive encounters with stars and 
Giant Molecular Clouds (GMCs) in the Galaxy, but are 
more likely imposed early on in the formation process. 

2. OBSERVATIONS 

We imaged a sample of ten T dwarfs identified in a 
magnitude-limited search of the 2MASS database (Bur¬ 
gasser et al. 1999, 2000b,c, 2002d) in the WPFC2 F814W 
(Ac = 7940A) and F1042M (Ac = 1020lA) filters during 
HST Cycle 9. A log of observations is given in Table 1. 
The F814W and F1042M filters were chosen to sample the 
red wing of the pressure-broadened K I doublet, as the 
strength of this feature allows late-type L and T dwarfs to 
stand out from nearly all background sources in red opti¬ 
cal colors (e.g., I-z; Leggett et al. 2000). Each object was 
centered on the PC chip and observed twice (to allow for 
cosmic-ray subtraction) in both filters for total exposure 
times ranging from 2000 to 2600 sec. 

Images were reduced by standard pipeline processing, 
which includes analog-to-digital correction, removal of the 
bias pedestal, subtraction of bias and dark frames, and 
division by an appropriate flat field image. No correction 
to shutter shading was required due to the long exposure 
times. The images were then combined using a cosmic 
ray rejection routine developed at NASA Goddard, based 
on an improved version of the CR_REJECT routine writ¬ 
ten by R. S. Hill. Bad pixels identified both in the data 


quality arrays and the cosmic ray rejection routine were 
replaced by the mean of neighboring pixels to produce the 
final, cleaned image. We note that a reflection arc from an 
offset bright star appears on the F814W PC chip images 
of 2MASS 2356—1553, caused by non-optimal baffling of 
this detector (Biretta et al. 2002). This reflection leads to 
a slightly higher background in the vicinity of the source 
but does not affect our background-subtracted photome¬ 
try (§3.1). No other residual or reflection artifacts are seen 
in the data. 

Sections of the reduced PC images 2'.'3 on a side 
around each of the primary targets are shown in Figure 1. 
North/east orientations are indicated by arrows. We im¬ 
mediately identify two closely-separated objects near the 
position of 2MASS 1225—2739 in both the F814W and 
F1042M images, with the western component appearing 
to be fainter at both bands. 2MASS 1534—2952 is slightly 
elongated along a north/south axis as compared to both 
the other T dwarf targets and other sources in the 2MASS 
1534—2952 PC field, as is shown in more detail in Figure 2. 
We attribute this elongation to a marginally resolved pair 
of point sources (§4.2). All of the other T dwarf targets 
appear to be single point sources at the spatial resolution 
of the PC chip (0"046). 

3. PHOTOMETRY 
3.1. Aperture Photometry 

Sources on all four WPFC2 chips were initially identified 
with the IRAF® DAOFIND routine, and then confirmed 
by visual inspection. We extracted aperture photometry 
for each source using the IRAF PHOT routine, using 2, 
3, and 5 pixel apertures, corresponding to (y.'09, 0'.'14, and 
0723 on the PC chip and 0720 , 0730, and 0750 on the WF 
chips. Background values were determined using a cen¬ 
troid algorithm in a 15 pixel annulus 10 pixels from each 
source. Individual aperture corrections were then mea¬ 
sured for each single star (neglecting galaxies and resid¬ 
ual cosmic rays) by integrating their radial profiles to 20 
pixels. Because our fields were generally at high galac¬ 
tic latitudes, source densities were low, and many fields 
(particularly those in the F1042M filter) had few or no 
stellar sources on a particular chip. Hence, a mean set of 
aperture corrections were derived from all point sources 
in each chip and filter; these values are given in Table 2. 
After applying the aperture corrections, flux values were 
corrected for geometric distortion and charge-transfer effi¬ 
ciency (GTE), the latter by the prescription of Whitmore, 
Heyer, & Casertano (1999). Because of the long exposure 
times and high backgrounds, typical GTE corrections gen¬ 
erally did not exceed 20%. No corrections for instrument 
contamination were made, as they are exceedingly small in 
these red filters (Biretta et al. 2002). Synthetic flux zero- 
points from Baggett et al. (1997) were used to convert the 
corrected magnitudes to the WPFG2 photometric system. 
A final source list was compiled by requiring detections in 
both filters and positional coincidence within 1 pixel; this 
constraint eliminated most residual cosmic rays (however, 
see §5.2). 


® Throughout the main text, we abbreviate object names to 2MASS hhmmiddmm; full designations are given in Table 1. 

® IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Because aperture corrections were made for point-source 
radial profiles, extended sources such as galaxies or close 
doubles were readily identified by significant differences in 
derived photometry depending on aperture size. This ef¬ 
fect verified the slight elongation of 2MASS 1534—2952. 
All other T dwarf targets (including the two close sources 
in the 2MASS 1225—2739 field) have photometry in each 
aperture consistent with the formal uncertainties. We 
adopt the 3-pixel aperture magnitudes for point sources 
(optimizing the signal-to-noise ratio), except for the sec¬ 
ond source in the 2MASS 1225—2739 field, where we se¬ 
lect a 2 -pixel aperture to minimize the contribution of the 
brighter nearby source; and extended sources, including 
2MASS 1534—2952, where we select 5-pixel aperture mag¬ 
nitudes to minimize aperture corrections. HST WPFC2 
and 2MASS J-band magnitudes and colors for our T dwarf 
targets and Gliese 229B^° (Golimowski et al. 1998; Leggett 
et al. 1999) are listed in Table 3. 

3.2. T Dwarf Colors 

Figure 3 plots F814W magnitude versus 
F814W—F1042M color for all sources identified in the 
ten WPFG2 datasets, along with data for Gliese 229B. 
Single point sources and target objects are plotted as solid 
circles, while extended sources (i.e., galaxies) are plotted 
as open circles. Primary T dwarf targets are individually 
labelled, all of which are 2-3 magnitudes redder than the 
background stellar and galactic sources, again due to the 
red wing of the pressure-broadened K I doublet. Both 
sources at the position of 2MASS 1225—2739 lie at red 
colors, implying that both are T dwarfs. Based on the es¬ 
timated surface density of T dwarfs detectable by 2MASS, 
8.4x10“^ deg“^ (Burgasser et al. 2002d), which we ex¬ 
trapolate to a limiting J magnitude of 17 (the apparent J 
magnitude of 2MASS 1225—2739B; see below), the prob¬ 
ability of two relatively bright T dwarfs randomly lying 
within ~18" of each other (the approximate search ra¬ 
dius on the PG chip) is 3x10“^. We therefore confidently 
claim companionship for these two objects based on their 
proximity and unique colors. By the same argument, the 
two sources at the position of 2MASS 1534—2352 are also 
companion T dwarfs, based on the red color of their com¬ 
bined light. Hereafter, we refer to these two systems as 
2MASS 1225-2739AB and 2MASS 1534-2352AB. 

Figure 4 plots the F814W—F1042M versus F814W—J 
color-color diagram for the observed T dwarfs and Gliese 
229B. Note that the colors of single targets follow a fairly 
linear trend: 

[J^8141F-J] = (3.1±0.3)-b(0.46±0.07)x [7^8141^-7^1042117]. 

( 1 ) 

Because 2MASS 1225—2739 is unresolved by 2MASS, we 
determined J-band component magnitudes from the com¬ 
bined light magnitude, J = 15.22±0.05, and the J-band 
flux ratio, 

AJ = AF8141F - A[7^8141F - J] = 1.35±0.08, (2) 

using Eqn. 1 and the photometry listed in Table 3. The 
F814W—J colors for these two objects both lie ~0.15 mag 
below the linear fit traced by the single stars, but are con¬ 
sistent within the photometric uncertainties. The com¬ 
bined light F814W-J color of 2MASS 1534-2952AB is 


also below the single star locus, but in this case it is 
probably because the object is marginally resolved in the 
WPFC2 images. On the other hand, 2MASS 1217—0311 
is slightly redder in F814W—J color than expected, al¬ 
though by no more than 2a. Both F814W—F1042M and 
F814W—J colors are generally redder for the later-type 
T dwarfs, with the former being particularly sensitive to 
spectral type. One notable exception is the T6.5 emission¬ 
line dwarf 2MASS 1237-1-6526 (Burgasser et al. 2000a), 
which has the reddest F814W—F1042M color in the sam¬ 
ple (see §5.3). On the other hand, F1042M—J colors gen¬ 
erally decrease for later spectral types, likely due to in¬ 
creased H 2 O and CH 4 absorption around 1.25 /im (Bur¬ 
gasser et al. 2002 d). 

4. PSF FITTING 

4.1. Technique 

In order to derive separations and flux ratios for our two 
T dwarf binaries, and search for faint companions around 
the other target sources, we performed point spread func¬ 
tion (PSF) subtraction on all of our primary targets. Our 
technique was as follows: first, we extracted subimages of 
all apparently single point sources from the PC chip im¬ 
ages of all ten datasets, a total of 22 sources in F814W 
and 11 in F1042M. These included some of the target ob¬ 
jects, although care was taken to exclude any point sources 
with bad pixels near the source peak. We then subtracted 
two-dimensional Gaussian fits to the PSFs from the im¬ 
ages; typical residuals were consistently < 10 % of the 
original source peak. Finally, we averaged these Gaussian- 
subtracted images, scaled by the fit maximum, to produce 
a single PSF residual image for each filter. 

For each of our target sources, we searched for faint com¬ 
panions using an iterative reduction routine. Model im¬ 
ages were constructed by combining two PSF residual im¬ 
ages with two Gaussian surfaces having the same FWHM 
as the PSF fits described above, each scaled to separate 
component fluxes. For 2MASS 1225—2739AB and 2MASS 
1534—2952AB, initial guesses were based on the approxi¬ 
mate positions and flux ratios from the aperture photom¬ 
etry (we assumed 2MASS 1534—2952AB to be separated 
by 1 pixel along each axis as an initial guess). Our routine 
then iteratively searched for the optimal solution to the 
primary position, secondary position, primary flux, and 
secondary flux, in that order, by shifting the component 
positions in steps of 0.1 pixels and scaling the fluxes in 
steps of 1% (0.01 mag). If the secondary flux was scaled 
below 1 count or separations below 0.5 pixels were reached, 
then the object was considered a single point source. Oth¬ 
erwise, the routine derived separations, position angles, 
and flux ratios for the optimal binary solution. 

For all of the apparently single targets, we followed up 
this process by fitting a single PSF residual plus Gaus¬ 
sian to the image and then searching by eye for any ob¬ 
vious counterparts. We then used the same binary search 
routine for each image with 20 random companion initial 
positions. If no companion brighter than the S/N = 7 de¬ 
tection limits (approximately 25.5 mag at F814W and 19.9 
mag at F1042M; Table 5) was found, the primary target 
was assumed to be single. 


J-band photometry for Gliese 229B is on the UKIRT system (Casali & Hawarden 1992), which is similar to the 2MASS photometric system 
(Carpenter 2001; Dahn et al. 2002). 
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4.2. Results 

Convergent binary solutions were obtained for both 
2MASS 1225-2739AB and 2MASS 1534-2952AB in both 
filters. Figure 5 shows the original and PSF-subtracted 
images for these pairs. Residuals from the subtraction 
were less than 4-5% of the peak flux, at the level of 8-10 
(F814W) and 2-3 (F1042M) times the background noise. 
In both cases, smaller residuals were obtained for fits to 
two sources rather than a single source. This validates 
the duplicity of 2MASS 1534—2952AB, whose separation 
(0'.'065±0'.'007) is smaller than both the diffraction limit 
(O'.'OS at F814W, O'.'ll at F1042M) and Nyquist sampling 
limit (2x0'.'046 = 0709) of the instrument. We are able 
to overcome the former because our technique resolves 
even significant overlap of two PSFs, particularly when 
they have nearly equal brightness. The latter constraint 
is overcome by using a PSF generated from multiple mea¬ 
surements, allowing us to subsample below the Nyquist 
limit. Table 4 lists the derived binary parameters for these 
systems; no other sources can be seen in the subtracted 
images. PSF fitting of the other primary T dwarf targets 
revealed only one potential faint companion to 2MASS 
1217—0311, detected at F1042M only. This possible de¬ 
tection, which may simply be a residual cosmic ray, is dis¬ 
cussed in detail in §5.2. No other companions were iden¬ 
tified around any of the other target objects within I'.'l, 
and no other faint objects with T dwarf-like colors were 
identified in any of the WPFC2 images. 

To obtain proper calibration and determine the uncer¬ 
tainties of our results, we ran the algorithm described 
above on 20,000 simulated binaries constructed from the 
F814W and F1042M images of 2MASS 0559-1404. These 
test images sampled a range of separations 1-15 pixels 
(0'.'05-0'.'69), all orientations, and flux ratios AM = 0- 
7 mag. Once processed through the PSF fitting routine, 
those test cases having output separations within 0.5 pixels 
and corrected flux ratios within 0.2 mag of the input values 
were considered recovered binaries. Corrections and uncer¬ 
tainties to both positions and flux ratios were then deter¬ 
mined separately for 2MASS 1225—2739AB and 2MASS 
1534—2952AB in each filter, using only those recovered 
test cases having similar input separations (a > 0723 and 
a < 0714, respectively) and flux ratios (1 < AM < 2.5 
and 0 < AM < 1, respectively). Typical flux ratio cor¬ 
rections were approximately — 0.10 mag (i.e., shifting the 
secondary to brighter magnitudes) with Icr uncertainties of 
0.04 and 0.3 mag for 2MASS 1225-2739AB and 2MASS 
1534—2952AB, respectively; separation Icr uncertainties 
were 0.12 and 0.15 pixels (07005 and 070 07), respectively, 
translating into position angle uncertainties of 7° and 9°. 
The values listed in Table 4 reflect these corrections and 
uncertainties. 

4.3. Search Limits 

Our calibration simulations allowed us to derive limit¬ 
ing detection magnitudes as a function of separation, as 
shown in Figure 6 . Around 2MASS 0559—1404, faint sec¬ 
ondaries (AM > 3) were generally missed at separations 
closer than 0715, while AM « 5.5 ( 6 ) and 4.5 (5) could be 
obtained for wide separations at the 95% (50%) confidence 


level at F814W and F1042M, respectively. In general, only 
near-equal-magnitude companions with a < 0709 could be 
recovered better than 50% of the time, as is the case for 
2MASS I534-2952AB. 

For a > Ol'd, S/N = 7 limits (Table 5) yield the max¬ 
imum sensitivity for faint companions, ranging from AM 
= 4.3-6.9 at F814W to AM = 2.9-4.9 at F1042M. We 
can convert these values to mass ratio {q = M 2 /M 1 ) lim¬ 
its using a mass-luminosity power law from Burrows et al. 
(2001), L oc and assuming for simplicity coevality 

and negligible variation in bolometric corrections over the 
sample^^, such that 

qiim = W ^ . (3) 

These values are listed in Table 5, and range from 1.0-0.4. 
Overall, our sample is complete to g > 0.4 for a > 4 AU 
(assuming a mean distance of 10 pc), with less sensitivity 
for small mass ratios to separations approaching a « 1 
AU. 

5. INDIVIDUAL TARGETS 
5.1. Binaries 

2MASS 1225—2739AB is clearly resolved into two 
unequal-magnitude components in our HST images. The 
colors of these objects are significantly different, with the 
fainter companion having an F814W—F1042M color simi¬ 
lar to the T8 Gliese 570D (Burgasser et al. 2000b), while 
the color of the primary is consistent with the spectral 
type of the combined system, T6. The magnitude ratios 
of this pair, AMf 8141 T = 1.59±0.04 and AMi?io 42 M = 
1.05±0.03, are greater than the absolute magnitude ra¬ 
tios of Gliese 229B (T6.5) and Gliese 570D, AMi? 8 i 4 VU = 
l.OliO.lO and AMi7’io42M = 0.87±0.12, consistent with 
2MASS 1225—2739A being earlier than type T6.5. Based 
on these colors, we speculate that this system is com¬ 
prised of a T6 and T8 pair, which should be confirmed 
with spatially-resolved spectroscopy. No parallax has been 
measured for this system yet, but the spectrophotomet- 
ric distance of the secondary, compared to Gliese 570D, 
is dfswiv = ^^fio 42 M = 11-1 pc- At J-band, dj = 10.8 
and dj = 11.9, based on the absolute J-band magnitudes 
of the T6 SDSS 1624-f0029 (Dahn et al. 2002, Mj = 
15.33±0.07) and Gliese 570D (Burgasser et al. 2000b, Mj 
= 16.47±0.05). All of these distance estimates combined 
yield a mean = 11.2±0.5 pc and projected physical 
separation a = 3.17±0.14 AU; note that the uncertainties 
do not represent probable scatter in the absolute magni¬ 
tude/spectral type relation, not yet adequately measured 
for the T dwarfs. Adopting Te// « 1000 and 800 K for 
the two components (Burgasser et al. 2002d), assuming 
coevality, and using the models of Burrows et al. (1997), 
we can derive component masses for ages of 0.5, 1.0, and 
5.0 Gyr, as listed in Table 6. The derived mass ratio of 
this system is q = 0.7-0.8, depending on its age. Assuming 
that, on average, the semimajor axis of a binary system 
< Usm >= 1.26 < a > (Fischer & Marcy 1992), we esti¬ 
mate orbital periods of 24-40 yr; hence, significant orbital 
motion (A(/) ~ 10°) could be detectable in this system on 
a yearly basis. 


The fact that F814W—F1042M varies by 1 mag over this sample implies that this is in general not true, particularly for very cool companions; 
however, this assumption is suitable for a rough limit. 
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2MASS 1534—2952AB is a more marginally resolved 
system, suggesting that it is much more closely-separated 
than 2MASS 1225—2739AB. The flux ratios for this sys¬ 
tem, AMp’ 8 i 4 n/ = 0.5i0.3 and AM 4 ?io 42 m = 0.2±0.3, 
are much smaller than the spread in absolute magni¬ 
tudes for mid-type T dwarfs; e.g., AMi? 8 i 4 iy = 2.06±0.06 
and AMfio 42 M = 1.36±0.08 between the T5 2MASS 
0559—1404 and the T6.5 Gliese 229B. Hence, we assume 
that this system is comprised of two nearly equal-mass 
components with spectral types T5.5 and Teff ~ HOO 
K. Again, no parallax measurement has been made for 
this system. Based on the absolute J-band magnitudes 
of 2MASS 0559-1404 and SDSS 1624-f0029, we estimate 
a spectrophotometric distance of dj = 16±5 pc, the un¬ 
certainty dominated by the poor absolute magnitude con¬ 
straints for mid-type T dwarfs. Combined with the mea¬ 
sured separation, this implies a = 1.0±0.3 A.U. Based on 
mass estimates as derived above (Table 6), we estimate or¬ 
bital periods of 4-6.5 yr, making this system an excellent 
target for rapid orbital mapping; however, its very close 
separation may hamper these measurements, and only sig¬ 
nificant orbital motion {A(f) > 30°) may be detectable. 

5.2. A Potential Binary 

PSF subtraction of the F1042M image of 2MASS 
1217—0311 reveal a faint companion 0'.'209±0'.'006 from 
the target source at position angle 74±7°. The corrected 
flux ratio of the secondary, AMi^’io 42 M = 2.35±0.04 makes 
this source the faintest “detection” in the sample, F1042M 
= 19.18±0.08, very close to the detection limits of the im¬ 
age. No counterpart is seen in the F814W image, suggest¬ 
ing a very red source, F814W-F1042M > 6.2. 

However, it is quite likely that this object is simply a 
residual cosmic ray. Figure 7 shows the two original PC 
exposures of the 2MASS 1217—0311 held, along with the 
hnal combined image, a single PSF-subtracted image, and 
a double PSF-subtracted image. Cosmic rays show up as 
bright pedestals of hux, in contrast to the more gradually 
sloping PSF. The potential companion (to the right of the 
primary target, indicated by an arrow in the fourth panel 
of Figure 7) is within one pixel of a bright, extended cosmic 
ray in the hrst image, and is completely wiped out by a cos¬ 
mic ray in the second image. It is likely that the source in 
the hrst image is a cosmic ray itself. Only 3-4 overlapping 
(i.e., detected on the same pixels in both exposures) cos¬ 
mic rays were found within 25 pixels of any of the target 
images, most of which were quite obviously identihable. 
This source is less obvious given its more gradually sloped 
prohle (i.e., not a hat pedestal) and very faint hux. Fit¬ 
ting two PSFs to the cosmic ray-corrected image results in 
signihcantly reduced residuals (by roughly 10%), although 
the fainter source is not completely subtracted out as a 
result of this ht (hfth panel in Figure 7). Such a ht is 
inconclusive, however, as subtracting any PSF from a cos¬ 
mic ray will reduce the overall image residuals, while faint 
legitimate sources may not subtract cleanly due to the in¬ 
creased relative noise. We therefore classify this source as 
a candidate companion, requiring additional followup to 
verify its existence. 

If the companion were real, it would be an extremely 
interesting object, as its color limit and intrinsic faintness 

Derivation of uncertainties for all sample statistics are described in 


would make it the coldest and faintest brown dwarf so 
far identihed. Assuming simplistically that logTg// scales 
with F1042M magnitude, and using the absolute magni¬ 
tudes and estimated Te//s of Gliese 229B (Marley et al. 
1996, ~950 K) and Gliese 570D (Burgasser et al. 2000b; 
Geballe et al. 2001, ~800 K), we estimate « 

0.6 and hence Pffj ~ 500 K. Because H 2 O begins to con¬ 
dense in atmospheres as cool as this (Burrows & Sharp 
1999), 2MASS 1217—0311B would probably not be a T 
dwarf but the prototype for a new spectral class. It would 
also have an extremely low mass, M « 0.012 Mq for an 
age of 1 Gyr (Burrows et al. 1997). Hence, conhrmation 
of this possible companion by follow-up imaging is clearly 
a priority. 

5.3. Single Sources 

Two other T dwarfs in our sample warrant additional 
discussion. The hrst is the bright (J = 13.83±0.03) T 
dwarf 2MASS 0559—1404 (Burgasser et al. 2000c). This 
object has a measured parallax (Dahn et al. 2002), and is 
over 1 mag brighter at J-band than the L8 dwarfs 2MASS 
1632-1-1904 (Kirkpatrick et al. 1999) and Gliese 584C 
(Kirkpatrick et al. 2000), and only 0.6 mag fainter at K^. 
Burgasser (2001) has suggested that 2MASS 0559—1404 
may be an equal-magnitude binary based on its bright¬ 
ness and derived Tg//, although an alternate hypothesis 
may be the rapid dissolution of dust cloud material over 
the L/T transition (Burgasser et al. 2002b). Our images 
rule out the presence of bright secondary closer than 0705. 
If this hypothetical companion exists, it is either currently 
aligned with the primary or is separated by less than 0.5 
AU. As at least one brown dwarf spectroscopic binary has 
been found with a separation less than this limit (Basri 
& Martin 1999, PPl 15), high-resolution spectroscopy of 
2MASS 0559—1404 may be required to fully rule out the 
presence of a close companion. 

2MASS 1237-1-6526 is another T dwarf whose duplicity 
is under consideration. Burgasser et al. (2000a) have sug¬ 
gested that the unique Ha emission in this object may be 
due to the presence of a close (a < 0.003 AU), interacting 
companion, although Burgasser et al. (2002a) have failed 
to Hnd photometric evidence of an eclipsing system. Our 
HST images do not rule out the presence of this hypotheti¬ 
cal companion, as a spatial resolution of 070002 (assuming 
a distance of 14 pc) would be required to resolve it. The 
very red F814W—F1042M and F814W—J colors of this 
T6.5 dwarf, similar to the T8 Gliese 570D, could arise from 
warm circum(sub)stellar material, consistent with this ob¬ 
ject being a very young and low-mass (3-12 Mj„p) weak- 
line T Tauri object (Burgasser et al. 2002a; Liebert et al. 
2002). Photometry from Dahn et al. (2002) conhrm this 
object’s very red optical/near-infrared colors, but give no 
evidence for a reddened J—K color (—0.26±0.20). Hence, 
the nature of this object remains ambiguous. 

6. BINARY STATISTICS 
6.1. The Binary Fraction 

Of the ten T dwarfs imaged in our sample, two have 
clearly resolved binary companions, implying an observed 
binary fraction of 20^7^%^^- This is consistent with re- 
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suits obtained for field late-type M and L dwarfs (Koerner 
et al. 1999; Reid et al. 2001; Close et al. 2002). 

One must keep in mind, however, that the T dwarf sam¬ 
ple was originally drawn from a magnitude-limited search, 
and that the observed binary fraction (for initially unre¬ 
solved pairs) is therefore biased (Opik 1924; Branch 1976). 
If we assume negligible contribution by multiple systems 
of three components or more, then the observed binary 
fraction, ^obs = is related to the “true”, or 

volume-limited, fraction, £& = Nunary/NtotaU by 


^obs 


1 - 66 - 1 - acb ’ 


where 


fo (1 + Pf^^fip)dp 
fo fip)dp 


( 4 ) 


( 5 ) 


is the fractional increase in volume sampled for binaries 
with flux ratio p = Fb/Fa and flux ratio distribution 
f{p). In the case of all binaries being equal-magnitude 
systems, a = 2 ^/^ = 2 . 8 ; while in the case of a flat /(p), a 
= 1.9. From these limiting cases and our observed binary 
fraction, we derive = 9 ^ 4 ^%, where we have included 
the uncertainty in and allowed a to vary between 1.9 
and 2.8. This value is significantly lower than the binary 
fraction of more massive stars (Duquennoy & Mayor 1991; 
Fischer & Marcy 1992; Reid & Gizis 1997). 

A second means of obtaining a bias-free estimate of 
the binary fraction in our magnitude-limited sample is 
by computing the relative luminosity functions with the 
l/Knaa; technique (Schmidt 1968, 1975). Simply, Cb = 
^binary/^total, where is $ is the luminosity function cal¬ 
culated from 


$ = 


N 


'^ = ^ y maa 


( 6 ) 


the sum carried over all N objects in the sample, with 
Vmax oc dmax ~ ^or the T dwarfs, the minimum 

detectable distance, dmin, is set by the constraint of no 
optical counterpart in the USNO A2.0 catalog (Monet et 
al. 1998) or in Digital Sky Survey images, and is roughly 
1 pc for all objects in our sample (Burgasser 2001); the 
maximum detectable distance, dmax, depends on the sam¬ 
ple search limit (J = 16) and the absolute J magni¬ 
tude of the object. We have estimated Mj for objects 
in our sample using the known absolute magnitudes of 
the T5 2MASS 0559-1404, T 6 SDSS 1624-^0029, T6.5 
Gliese 229B, and T 8 Gliese 570D, interpolating by spec¬ 
tral type (Burgasser 2001). Because binaries can be de¬ 
tected to distances yT^Fp further than single objects, we 
have included this correction for 2MASS 1225—2739AB 
and 2MASS 1534-2952AB using pj = 0.29 and 1.0, re¬ 
spectively. The derived binary fraction is only 6 %, on the 
low end of, but not inconsistent with, the bias-corrected 
value given above, and again much lower than stellar bi¬ 
nary fraction measures. An estimate of uncertainty for this 
technique is not straightforward (Mendez & Ruiz 2001); 
nonetheless, the value is consistent with a binary fraction 
much less than that of more massive stellar systems. The 
completeness estimator for our sample (Schmidt 1968), 
< V/Vmax > = 0.51±0.09, gives us some confidence that 
our result is not significantly influenced by incompleteness 
or color bias. 


The bias correction given in Eqn. 4 is applicable to the L 
dwarf sample of Reid et al. (2001), which is also based on 
a magnitude-limited survey of the 2MASS database (Kirk¬ 
patrick et al. 1999, 2000). They found ej)**® = 20^6^% (4 of 
20 ), which translates into a corrected fraction et, = 9 ^ 4 ^%, 
where again we have included uncertainty in both e^**® and 
a. The llVmax technique gives a consistent value of 12%; 
again, significantly lower than the binary fraction of more 
massive stars. 

Hence, we find that both L and T dwarf samples, when 
corrected for selection bias, yield binary fractions which 
are significantly lower than measurements made for more 
massive stars, suggesting a continuation of the apparent 
trend of decreasing et from F-G to M dwarf stars. How¬ 
ever, it must be stressed that the derived fraction is appli¬ 
cable for separations a > 1 — 5 AU and q > 0.4, while 
the investigation of, e.g., Duquennoy & Mayor (1991), 
probed much smaller mass ratios (q 0 . 1 ) and separa¬ 
tions (a ^ 0.1 AU). Hence, our sample may contain bi¬ 
naries with secondaries below our detection limits, or very 
tight unresolved binaries. We can estimate the contribu¬ 
tion from the latter population by examining the frequency 
of M dwarf spectroscopic binaries in the magnitude-limited 
sample of Reid et al. (2002), who found 6 ^ 2 % (2 of 36) 
of their targets were spectroscopic binaries. Again, us¬ 
ing Eqn. 4, this implies a bias-corrected fraction of only 
3 ^ 2 %) increasing the net binary fraction of the T dwarfs 
to perhaps 12 %, not enough to bring our results in agree¬ 
ment with the binary fraction of E-G or M dwarfs. Simi¬ 
larly, if we compare our derived binary frequency to only 
those stellar systems having a > 5 AU, roughly 41% in 
the Duquennoy & Mayor (1991) F-G dwarf sample (using 
their Gaussian log(P) distribution) and 31% in the Fis¬ 
cher & Marcy (1992) M dwarf sample (their Table 2), we 
find that there are clearly fewer T dwarf multiple systems 
in this separation regime. The contribution of lower-mass 
companions is also insufficient to explain the deficiency of 
L and T dwarf binaries, as binary fractions for F-G and M 
dwarf systems with 0.4 < q < 1.0 are roughly 33% (from 
Table 7 in Duquennoy and Mayor 1991) and 32% (from 
Fischer & Marcy 1992, assuming a flat mass ratio distri¬ 
bution), respectively, significantly higher than our results. 
Hence, unless T dwarfs prefer very closely-separated (see 
below) and/or very low-mass companions, the binary frac¬ 
tion of these objects is significantly lower than that of more 
massive stars. 


6.2. Separation Distribution 

The two confirmed binary systems identified in this sur¬ 
vey have projected separations a < 3 AU, and no wide, co¬ 
moving companions to any of these objects have yet been 
identified in either the HST data or the 2MASS survey. 
In fact, no wider companions (a > 2") have been found 
around any T dwarf identified in the 2MASS or SDSS sur¬ 
veys. This result is consistent with the current absence of 
widely-separated late-M and L dwarf binaries (Table 7), 
all of which have a < 10-20 AU. In contrast, roughly 50% 
of the more massive early-type M dwarf multiple systems 
in the Fischer & Marcy (1992) study have 10 AU a < 
10"^ AU. Similarly, roughly 40% of M dwarf multiple sys¬ 
tems in the 8 pc sample have separations greater than 10 
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AU (Reid & Gizis 1997). If lower-mass systems had a sim¬ 
ilar fraction of wide binaries, then roughly 20 pairs with 
a > 10 — 20 AU from the approximately 300 known L and 
T dwarfs should have been identified, while there are cur¬ 
rently none. The absence of wide systems may contribute 
to the overall dehciency in multiple systems amongst the 
T dwarfs, as the binary fractions of F-G dwarfs and M 
dwarfs drop to roughly 20% for separations a < 10 AU, 
within the uncertainty estimates of our bias-corrected re¬ 
sult. We discuss the apparent limit in the separations of 
low-mass stars and brown dwarfs further in §7.1. 

6.3. Mass Ratio Distribution 


trinsic faintness of very low-mass brown dwarfs, and the 
preferential selection of equal-mass systems in magnitude- 
limited surveys. We must, therefore, consider selection 
biases in these magnitude-limited samples. The correction 
to the flux ratio distribution is 


fip) 


(X (1-bp)^/^ 


(7) 


which, using the mass/flux ratio scaling as before, yields 


/(<?) 


Oc(l + g2.64)3/2^ 


( 8 ) 


Finally, we consider the mass ratio distribution, f{q), a 
statistic that can constrain the origin of secondaries in a 
binary population. In general, masses are difficult to de¬ 
rive for field brown dwarfs, as estimates depend on both 
temperature and age, and there are few empirical clues 
currently known for the latter parameter. In Table 7, 
we have estimated masses for 2MASS 0746-I-2000AB, DE¬ 
NIS 1228-II59AB, 2MASS 1534-2952AB, and 2MASS 
1225—2739AB assuming an age of 1 Gyr, the Tg// scale 
of Burgasser et al. (2002d), and the theoretical models of 
Burrows et al. (1997); for 2MASS 1146-I-2230AB, we used 
maximum masses of 0.06 M© based on the presence of the 
6708 A Li I line in the combined light spectrum (Kirk¬ 
patrick et al. 1999). All other mass estimates are taken 
from the listed references. Fortunately, the desired quan¬ 
tity, q, is not greatly sensitive to these assumptions. 

The two T dwarf binaries identihed in our sample have 
relatively large mass ratios, q = 0.8 and 1.0. As discussed 
in §4.3, we were capable of identifying systems down to 
q = 0.4, albeit not for very closely separated systems like 
2MASS 1534—2952AB. When we place these two systems 
in context with the other low-mass binaries listed in Table 
7, there appears to be a preference for equal-mass systems, 
the lowest mass ratio system being 0.7. This is similar 
to what has been reported in the 8 pc sample (Reid & 
Gizis 1997), and is at odds with the flatter distributions of 
Duquennoy & Mayor (1991) and Fischer & Marcy (1992). 

The two binaries identified in our program form too 
small a sample to examine the mass ratio distribution 
statistically, so we combined our ten T dwarfs with the 
L dwarf sample of Reid et al. (2001). Based on their 
completeness limits, and using the same mass/flux ratio 
scaling as described in §4.3, we find that their sample is 
complete to <7 > 0.4 for a > 0'.'25, or a > 5 AU assum¬ 
ing a mean distance of 20 pc. This is comparable to our 
completeness for a > 4 AU, although the inner separa¬ 
tion limit for our sample is roughly one-half that of the L 
dwarf study. Nonetheless, because similar instruments and 
observing strategies were employed, combining these two 
samples should not introduce signihcant biases. The ob¬ 
served binary fraction for this combined sample is 20 ^®%, 
while the bias-corrected fraction is 9^®%. Breaking these 
systems down by mass ratio, we find systems with 

1.0 < g < 0.9, 11Lq 4 system with 0.9 < <7 < 0.8, I1Lq 4 
system with 0.8 < q < 0.7, and less than 1.4 systems for 
all other ratio bins. We plot this distribution (light grey 
histogram) in Figure 8 , normalized so that f{q = 1) = 1. 

Again, there appears to be a preference for equal mass 
binaries. This is not unexpected, however, given the in¬ 


Hence, the bias is a fairly strong function of q. Apply¬ 
ing these corrections, we derive a slightly revised mass 
ratio distribution (dark grey histogram in Figure 8 ). Even 
with the bias corrections, there are more objects in the 
0.9 < q < 1.0 bin than in other mass ratio bins, although 
not a statistically signihcant number. With the substantial 
statistical uncertainties of our small sample, in particular 
the large upper limits for q < 0 . 6 , we cannot statistically 
rule out a flatter distribution. 

We have also plotted the M dwarf f{q) for the 8 pc sam¬ 
ple (Reid & Gizis 1997) in Eigure 8 , which we have nor¬ 
malized as above and computed uncertainties based on a 
total of 21 M dwarf multiple systems. This distribution ap¬ 
pears to be similar to the L and T dwarf distribution, with 
a possible preference for high-mass ratio systems, although 
the uncertainties are again signihcant; a hatter distribu¬ 
tion cannot be statistically ruled out. Hence, concluding 
a preference of equal-mass components amongst the M, L, 
and T dwarf binaries requires considerably better statis¬ 
tics, but our results are suggestive of this trend. 

7. BROWN DWARF BINARY FORMATION AND DISRUPTION 

The results above indicate that both the binary frac¬ 
tion and separation distribution of brown dwarfs are sig- 
nihcantly diherent that those of more massive stars, while 
the mass ratio distribution suggests a preference for equal- 
mass systems. We now examine how these properties may 
constrain the formation or evolution of substellar binary 
systems. 

7.1. Disruption by Stellar and GMC Encounters 

The dehciency of brown dwarf binaries with a > 10 AU 
is reminiscent of the dehciency of stellar binaries with a 
> 0.1 pc « 2x10® AU (Bahcall & Soniera 1981; Glose, 
Richer, & Crabtree 1990; Wasserman & Weinberg 1991). 
While there remains some debate as to whether a sharp 
break exists in the separation distribution (Retterer & 
King 1982; Wasserman & Weinberg 1987, 1991; Close, 
Richer, & Crabtree 1990), it is generally believed that im¬ 
pulsive perturbations by close stellar encounters or passage 
through a GMC causes a gradual diffusion of separations 
and binding energies, ultimately resulting in the dissolu¬ 
tion of weakly bound systems in a catastrophic encounter 
(Weinberg, Shapiro, & Wasserman 1987). Because the 
binding energies of brown dwarf pairs are small and such 
systems therefore easily disrupted, it is tempting to ascribe 
the same mechanism to the apparent absence of widely- 
separated systems. 
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To examine the probability of disruption by stellar 
and GMC encounters, we used the formalism of Wein¬ 
berg, Shapiro, & Wasserman (1987), adopting the gen¬ 
eral parameters used by the authors^^ and examine the 
evolution of two 0.05 Mq gravitationally bound brown 
dwarfs with separation 10 AU. The critical impact pa¬ 
rameter for significant gravitational disruptive effects is 
bmax oc « 70 AU. For stellar encounters, 

the Fokker-Planck impact parameter in the tidal limit^^, 
b*pp oc K, 30 AU « a, implies that both 

close, catastrophic collisions and gradual tidal disruption 
can affect the evolution of brown dwarf pairs. How¬ 
ever, the frequency of close stellar encounters is F*^; oc 
aM~^ « (2x10® Gyr)“^, while the diffusive timescale is 
T* oc a~^M « 3700 Gyr. Hence, stellar encounters are 
not frequent enough to affect brown dwarf binaries with 
a < 10'^ AU over the age of the Galaxy. The tidal limit 
impact parameter for GMC interactions is « 2x10® 

AU >> bmax, while the impact parameter for catas¬ 
trophic interactions with GMC clumps in the case of cloud 
penetration (occurring at a rate of roughly 1 Gyr“®) is 
^ 2000 AU >> bmax', hence, GMC interactions 
play no role in the disruption of brown dwarf binary sys¬ 
tems. Therefore, the separation limit of brown dwarf pairs 
is not due to disruption in the Galactic field, as appears to 
be the case for wide stellar binaries. Indeed, only brown 
dwarf systems with separations many orders of magnitude 
larger than those observed could be disrupted in the field. 

To further elucidate how stellar and GMC disruptions 
do not constrain the separation of brown dwarf pairs, Fig¬ 
ure 9 plots the separation of binary stars and brown dwarfs 
versus total mass. Binary data for brown dwarf and late- 
type stars (primaries later than M8) are listed in Table 7; 
for stellar binaries, we include the samples of Close, Richer, 
& Crabtree (1990), Duquennoy & Mayor (1991), Fischer & 
Marcy (1992), Reid & Gizis (1997, HST M dwarf binaries 
and the 8 pc sample), and Tokovinin (1997, Multiple Star 
Catalog); finally, for stellar-brown dwarf binaries we use 
compiled values from Reid et al. (2001). The absence of 
wide low-mass pairs is quite striking in this figure, particu¬ 
larly given the ability of 2MASS and other field surveys to 
resolve such systems. The curved line shows a log-normal 
relation for the maximum separation of binary systems, 
logOmaa: = 3.33Mtot +1.1, which is appropriate for disrup¬ 
tion by point-source encounters (Reid et al. 2001). Note 
that such an envelope matches the observed cut-off for stel¬ 
lar binaries quite well, but allows more-widely separated 
brown dwarf binaries to form (a « 20 — 30 AU). For the 
lower-mass systems, we find a second line, 

amaxiAU) = UOOxMIp (9) 

appears to be more adequate for the separation limit. 
While the number of objects for which this envelope ap¬ 
plies is relatively small, it is not biased by selection effects, 
as all systems would be unresolved in their original sur¬ 
veys, and would be easily resolved by HST, for a > Gmax- 
We suggest that this power-law relation may be a clue to 
the mechanism that modulates the formation or disruption 
of substellar binaries, although further data are required to 
confirm if this relation is truly representative of all brown 


dwarf binary systems. 

What about disruption within the nascent star-forming 
cluster? A survey of the ~ 120 Myr Pleiades cluster by 
Martin et al. (2000) found no binaries out of a sample 
of 34 low-mass star or brown dwarf members for a > 27 
AU, although candidate photometric binaries (including 
PPL 15) suggest a binary fraction of ~ 22% (Reid et al. 
2001 ), consistent with the binary fraction observed in the 
field. Duchene, Bouvier, & Simon (1999) found no sub- 
stellar binaries or companions in the 0.5-10 Myr IG 348 
cluster for a > 30 AU, although one very wide (a « 2300 
AU) candidate system in this cluster has been suggested 
by Najita, Tiede, & Carr (2000). Finally, a search for bi¬ 
nary objects in 1-5 Myr Cha Ha 1 cloud by Neuhauser et 
al. (2002) and Neuhauser, Guenther, & Brandner (2003) 
has turned up only one potential binary candidate with 
a < 28 AU. Therefore, it appears that for disruption to 
play an important role in the elimination of brown dwarf 
binaries with a > 10 AU, it must occur within a few mil¬ 
lion years of formation. Note that the theory of Wein¬ 
berg, Shapiro, & Wasserman (1987) predicts that stellar 
encounters may have some influence in young dense clus¬ 
ters, as F*^( oc n* ~ 1 Gyr“^, and r* oc ~ 20 Myr for 
n ~ 10^ pc“®, typical for regions such as the Orion Neb¬ 
ular Cloud (Hillenbrand 1997). However, as these dense 
regions rapidly disperse (i.e., within a few Myr), close en¬ 
counters are probably not solely responsible for the ab¬ 
sence of widely-separated substellar binaries. 


7.2. Small N Protoclusters and Brown Dwarf Ejection 

Models 

A currently popular model of star formation in clusters 
is through the fragmentation of molecular clouds into small 
aggregates of non-hierarchical protostellar cores (Larson 
1972), with the entire young star forming region being 
comprised of these initial groupings. On a short timescale 
(ages < 10® yr), these “protoclusters” are disrupted by 
dynamical interactions between the cores, which is likely 
modulated by residual gas and dust that continues to be 
accreted (Bonnell et al. 1997, 2001). Such dynamical in¬ 
teractions preferentially eject the lowest-mass components, 
while an ejected core is also less likely to continue sig¬ 
nificant accretion. These considerations have given rise 
to so-called “ejection” formation models for brown dwarfs 
(Reipurth & Glarke 2001; Bate, Bonnell, & Bromm 2002), 
in which the dynamic removal of cores from accretion re¬ 
gions condemns them to remain below the Hydrogen burn¬ 
ing minimum mass. A numerical simulation by Bate, Bon¬ 
nell, & Bromm (2002) utilizing this general model have 
found a low brown dwarf binary fraction, at most 5%, 
based on a single remaining undisturbed pair in a dynam¬ 
ically unstable multiple system. This fraction is consis¬ 
tent with the derived fraction of L and T dwarf bina¬ 
ries, although the simulation also predicts similarly low 
binary fractions for low-mass stars, which is not observed. 
Nonetheless, since the ejection model predicts the disrup¬ 
tion of potential brown dwarf binaries at very early ages, 
while also imposing a limit to the dimensions of such sys- 


V^ei = 20 km s e = 0.1, n* = 0.05 pc riQMC = 4x10 ® pc Rgmc = 20 pc, Mqmc = 5x10® Mq, and = 25; see Weinberg, 

Shapiro, & Wasserman (1987) for nomenclature. 

In the case we are considering, GM/eaVf^i > (M/M*) ~ 0.1, so that, unlike wide stellar pairs, the tidal limit applies. 
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terns (Reipurth & Clarke 2001), it shows some promise in 
explaining the origins of substellar systems in general. 

7.3. Fragmentation 

The preference for brown dwarfs to form close binaries 
may not necessarily require a disruptive process, however. 
Studies of young binary stars favor fragmentation (Boss 
1988) as the dominant mode of binary formation, due 
to coevality of components, the presence of circumbinary 
structures, and the preference for equal-mass components 
in closely-separated systems (White & Ghez 2001). These 
conditions do not require dynamical disruption from neigh¬ 
boring protostellar systems. In general, a low-mass gas 
and dust core must collapse to smaller dimensions before 
it achieves sufficient densities to continue fragmentation, 
producing multiple systems which are initially closely sep¬ 
arated. This suggests a maximum separation dependence 
on mass, as hinted at in Figure 8, although no theoreti¬ 
cal prediction as such has been made. The deficiency of 
low-mass pairs may arise from the inability for very small 
cloud clumps to both form and also continue fragmenting, 
although the influence of magnetic fields, turbulence, and 
external perturbations would also have substantial influ¬ 
ence. Current models (e.g.. Boss 2001) are capable of pro¬ 
ducing core fragments in the range of 10s of Jupiter masses 
(Mjup), in the mass range of field L and T brown dwarfs, 
but masses down to 1 Mjup require dynamical ejection to 
prevent further accretion. 

7.4. How Do Brown Dwarfs Form? 

The similarity in the binary fractions and separation 
distributions for young cluster and field low-mass systems, 
and the low probability of dynamic disruption in all but 
the densest stellar environments, makes it highly proba¬ 
ble that the field brown dwarf binary distribution is quite 
similar to the natal distribution. This is important, as the 
distances and dust opacity of protostellar environments, 
and the relative faintness of protosubstellar objects, makes 


investigation of brown dwarf formation at very early ages 
quite difficult. Improving the statistics for field brown 
dwarf systems, and examining closer separation regimes 
through radial velocity techniques, should provide consid¬ 
erable insight into the formation of these very low-mass 
objects. We find that both fragmentation and ejection 
models produce some of the qualitative characteristics of 
late-M, L, and T dwarf binaries, and it is possible that 
substellar systems form by some combination of these pro¬ 
cesses. However, more detailed quantitative predictions 
must be matched with large, unbiased sample statistics 
before conclusive statements can be made on the forma¬ 
tion of brown dwarfs. 
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APPENDIX 


PROBABILITY DISTRIBUTION FOR THE BINARY FRACTION 

When binary fractions (or other equivalent frequency statistics) are quoted in the literature, they are frequently assigned 
Poisson uncertainties. However, the Poisson limit applies only in the case of a large sample, whereas the brown dwarf 
samples discussed here are less than 30 in number. Hence, we derived statistical uncertainties by constructing a probability 
distribution for Cf, given the total sample size, N, and the number of binaries in the sample, n. The binomial distribution 
determines the probability of finding n binaries given the sample size and binary fraction, as: 

However, this equation may also be used to derive the probability distribution of es given the observed quantities N and 
n. To do this, we compute Bfei,; n, N) cx B{n; N, et) for 0 < Cf, < 1, normalizing 


f B'{eb-,n,N)deb = I, 

JO 


(A2) 


which yields B' = {N + 1)B. 

Figure 9 plots B' for our T dwarf sample, = 10 and n = 2. To derive upper and lower uncertainty limits, and 
e^, we computed the values for which B'deb = f^L B'deb = 0.84, equivalent to Icr limits for a Gaussian distribution. 
These limits can also be found numerically by solving 

(N+iy. 0M,x = ei 

0.16, X = 


E 

i=0 


d(Af-bl-f)! 


x\l-xy 


(A3) 
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As shown in Figure 9, the derived limits are not symmetric about the probability peak, prohibiting ranges which exceed 
the sample size or are less than zero. For large samples {N > 100) one recovers the standard Poisson uncertainty limits, 
- eb)/eb = (e& - e^)/e& = i/l/n + 1/N. 
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Table A1 

Log of HST Observations. 


Object^ 

(1) 

SpT 

(2) 

F814W 


F1042M 


RA (°)° 
(7) 

UT Date/Time"^ 
(3) 

t (sec) 

(4) 

UT Date/Time"^ 
(5) 

t (sec) 

(6) 

2MASS 05591914-1404488 

T5 

20000906 18:41 

2400 

20000906 20:15 

2600 

305 

2MASS 09373487+2931409 

T6pec 

20001016 23:44 

2600 

20001017 00:08 

2600 

339 

2MASS 10475385+2124234 

T6.5 

20010104 16:16 

2400 

20010104 17:50 

2600 

310 

2MASS 12171110-0311131 

T7.5 

20000704 03:10 

2400 

20000704 04:43 

2600 

158 

2MASS 12255432-2739466 

T6 

20010410 15:53 

2600 

20010410 17:29 

2600 

79 

2MASS 12373919+6526148 

T6.5 

20000613 15:32 

2000 

20000613 16:49 

2400 

160 

Gliese 570D 

T8 

20000818 07:23 

2400 

20000818 08:54 

2600 

152 

2MASS 15344984-2952274 

T5.5 

20000818 04:10 

2600 

20000818 05:42 

2600 

148 

2MASS 15462718-3325111 

T5.5 

20000819 05:55 

2600 

20000818 07:26 

2600 

147 

2MASS 23565477-1553111 

T6 

20001129 09:48 

2400 

20001129 11:23 

2600 

66 


^'Source designations for the 2MASS Point Source Catalog are given as “2MASS Jhhmmss[.]ss±ddmmss[.]s”. The 
suffix conforms to lAU nomenclature convention and is the sexigesimal R.A. and decl. at J2000 equinox. 

^UT date/time given as yyyymmdd hh:mm. 

‘^Telescope Roll Angle East from North. 


Table A2 

WPFC2 Aperture Corrections. 


Filter 

(1) 

Aperture^ 

(2) 

PC 

(3) 

WFl 

(4) 

WF2 

(5) 

WF3 

(6) 



24 

66 

79 

78 


2 

1.62+0.07 

1.21+0.03 

1.21+0.02 

1.21+0.02 

F814W 

3 

1.21+0.04 

1.09+0.02 

1.08+0.01 

1.08+0.01 


5 

1.04+0.02 

1.02+0.02 

1.01+0.01 

1.01+0.01 


N* 

7 

21 

17 

19 


2 

1.65+0.08 

1.28+0.06 

1.31+0.05 

1.28+0.05 

F1042M 

3 

1.33+0.09 

1.11+0.04 

1.13+0.03 

1.12+0.03 


5 

1.01+0.06 

1.03+0.03 

1.03+0.02 

1.03+0.02 


Aperture in pixels, corresponding to angular apertures of 0^^09 (0(^20), 0(T4 
(0^/30), and 0^.^23 (0^.^50) for 2, 3, and 5 pixels on the PC (WF) chips. 


Table A3 

T Dwarf Photometry. 


Object 

(1) 

SpT 

(2) 

F814W 

(3) 

F1042M 

(4) 

2MASS J 
(5) 

F814W-F1042M 

(6) 

F1042M-J 

(7) 

F814W-J 

(8) 

2MASS 0559-1404 

T5 

18.65+0.03 

15.02+0.07 

13.83+0.03 

3.64+0.08 

1.19+0.08 

4.82+0.04 

2MASS 1534-2952AB'* 

T5.5 

19.62+0.02 

15.75+0.07 

14.90+0.04 

3.86+0.08 

0.85+0.08 

4.72±0.04 

2MASS 1546-3325 

T5.5 

20.52+0.03 

16.66+0.07 

15.60+0.05 

3.86+0.08 

1.06+0.09 

4.92+0.06 

2MASS 1225-2739A'’ 

T6 

20.32+0.03 

16.38+0.07 

15.50+0.05 

3.94+0.08 

0.88+0.09 

4.83+0.06 

2MASS 2356-1553 

T6 

20.73+0.03 

16.96+0.08 

15.80+0.06 

3.77+0.08 

1.16+0.10 

4.93+0.07 

2MASS 0937+2931 

T6pec 

19.73+0.03 

15.47+0.07 

14.65+0.04 

4.26+0.08 

0.82+0.08 

5.08+0.05 

Gliese 2293*^ 

T6.5 

19.49+0.03 

15.16+0.03 

14.32+0.05 

4.33+0.04 

1.15+0.06 

5.17+0.06 

2MASS 1047+2124 

T6.5 

20.89+0.03 

16.67+0.07 

15.82+0.06 

4.22+0.08 

0.85+0.10 

5.07+0.07 

2MASS 1237+6526 

T6.5 

21.25+0.03 

16.76+0.07 

16.03+0.09 

4.49+0.08 

0.73+0.12 

5.22+0.10 

2MASS 1217-0311 

T7.5 

21.11+0.03 

16.83+0.07 

15.85+0.07 

4.29+0.08 

0.98+0.10 

5.27+0.08 

Gliese 570D 

T8 

20.55+0.03 

16.08+0.07 

15.33+0.05 

4.49+0.08 

0.73+0.09 

5.22+0.06 

2MASS 1225-2739B*’''^ 

T8:® 

21.91+0.05 

17.45+0.06 

16.85+0.08 

4.46+0.07 

0.61+0.10 

5.07+0.09 


®^WPFC2 magnitudes computed for a 5-pixel aperture. 

*^J magnitude estimated from combined 2MASS J — 15.22±0.05 and AJ — 1.35±0.08 (see §3.2). 
‘^WPFC2 magnitudes from Golimowski et al. (1998); UKIRT J from Leggett et al. (1999). 
‘^WPFC2 magnitudes computed for a 2-pixel aperture. 

^Spectral type estimated from F814W—F1042M color; see §5.1.1. 
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Table A4 

Binary Parameters. 



2MASS 1225-2739AB 
(1) 

2MASS 1534-2952AB 
(2) 

2MASS 1217-0311AB?'‘ 
(3) 

SpT 

T6/T8 

T5.5/T5.5 

T7.5/Y? 

d (pc)'’ 

11.2±0.5 

16±5 

10±4 

a (") 

0"282±0"005 

0"065±0"007 

0'.'209±0'.'006 

a (AU) 

3.17±0.14 

1.0T0.3 

2.1±0.8 

<!> n 

250±7 

1±9 

74±7 

AF814W 

1.59±0.04 

0.5±0.3 

> 4.4 

AF1042M 

1.05±0.03 

0.2±0.3 

2.35±0.04 


^Potential faint companion requiring confirmation; see §5.1.2. 

^Spectrophotometric distance estimated from spectral types and T dwarf with known dis¬ 
tances; see §5.1.1. 


Table A5 

Limiting Detection Magnitudes for a > 0'.'4. 


Object 

(1) 


F814W 



F1042M 


^lim 

(2) 

AMlim 

(3) 

Qlim 

(4) 

(5) 

AMlim 

(6) 

Qlim 

(7) 

2MASS 0559-1404 

25.6 

6.9 

0.09 

19.9 

4.9 

0.18 

2MASS 0937+2931 

25.4 

5.8 

0.13 

19.9 

4.4 

0.22 

2MASS 1047+2124 

25.5 

4.6 

0.20 

19.9 

3.2 

0.33 

2MASS 1217-0311 

25.6 

4.4 

0.22 

19.9 

3.1 

0.34 

2MASS 1225-2739 

25.5 

5.2 

0.16 

20.0 

3.5 

0.29 

2MASS 1237+6526 

25.3 

4.1 

0.24 

19.7 

3.1 

0.34 

Gliese 570D 

25.4 

5.0 

0.17 

20.0 

3.8 

0.27 

2MASS 1534-2952 

25.2 

5.9 

0.13 

19.6 

4.2 

0.23 

2MASS 1546-3325 

25.4 

5.0 

0.17 

19.9 

3.2 

0.33 

2MASS 2356-1553 

25.4 

4.8 

0.19 

20.0 

2.9 

0.36 


^S/N = 7 detection limit. 

^Mass ratio limit derived from Eqn. 3. 


Table A6 

Estimated Orbital Parameters. 


Age (Gyr) 
(1) 

2MASS 1225-2739AB 

2MASS 1534-2952AB 

M (Mo) 

(2) 

a,m (AU)'‘ 
(3) 

P (y) 

(4) 

M (Mo) 

(5) 

a,m (AU)“ 
(6) 

P (yr) 

(7) 

0.5 

0.017/0.023 

4.0 

40 

0.027/0.027 

1.3 

6.4 

1.0 

0.024/0.033 

4.0 

34 

0.035/0.035 

1.3 

5.6 

5.0 

0.04/0.06 

4.0 

24 

0.065/0.065 

1.3 

4.1 


Semimajor axis assuming < asm > = 1.26< a > (Fischer Sz Marcy 1992). 
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Table A7 

Field Late-M, L, and T Dwarf Binaries. 


Object 

(1) 

SpT 

(2) 

M (Mq) 

(3) 

q 

(4) 

a (") 

(5) 

a (AU)“ 

(6) 

AM 

(7) 

Filter 

(8) 

Ref. 

(9) 

2MASS 2206-2047AB 

M8/M8.5 

0.090/0.088 

1.0 

0"17 

4.1 

0.11 

JHK^K 

1 

Gliese 569BC 

M8.5/M9 

0.069/0.059 

0.9 

0('08 

0.90±0.02 

0.45 

JK 

2 

2MASS 2331-0406AB 

M8/L3 

0.091/0.062 

0.7 

0"57 

14.4 

2.6 

JHK'K 

1 

LHS 2397aAB 

M8/L7.5 

0.090/0.068 

0.8 

0''21 

3.0 

4.5 

F814W 

3 

2MASS 1426+1557AB 

M8.5/L1 

0.083/0.075 

0.9 

0('l5 

3.6 

0.67 

JHK^K 

1 

2MASS 2140+1652AB 

M8.5/L0 

0.087/0.075 

0.9 

0('l6 

3.7 

0.75 

JHK^K 

1 

2MASS 0746+2000AB 

L0.5/L0.5 

0.075/0.075*’ 

1.0 

0('22 

2.7 

0.63 

F814W 

4 

2MASS 1146+2230AB 

L3/L3 

0.06/0.06° 

1.0 

0(^29 

7.6 

0.31 

F814W 

4 





0(^29 

7.6 

0.0 

K 

5 

Gliese 564BC 

L4/L4 

0.053/0.053 

1.0 

0"13 

2.4 

0.30 

JHKs 

6,7 

DENIS 1228-1159AB 

L5/L5 

0.06/0.06° 

1.0 

0^^28 

5.1 

0.22 

FllOM 

8,4 





0"27 

4.9 

0.10 

K 

5 

2MASS 0850+1057AB 

L6/T? 

0.05/0.04 

0.8 

o'.'w 

4.4 

1.3 

F814W 

4 

2MASS 0920+3517AB 

L6.5/L6.5 

0.68/0.68 

1.0 

0'.'07 

1.6 

0.43 

F814W 

4 

DENIS 0205-1159AB 

L7/L7 

0.05/0.05 

1.0 

O'.'dl 

9.2 

0.0 

K 

5 





0'.'35 

6.3 

0.0 

JHKL’ 

9 

2MASS 1534-2952AB 

T5.5/T5.5 

0.035/0.035*’ 

1.0 

0"07 

1.0 

0.5 

F814W 

10 

2MASS 1225-2739AB 

T6/T8 

0.033/0.024*’ 

0.7 

0^^28 

3.2 

1.6 

F814W 

10 


^'Projected separation, except for Gliese 569BC whose orbit has been mapped (Lane et al. 2001). 

^Assuming and age of 1 Gyr, Te// scale from Burgasser et al. (2002d), and evolutionary models from Burrows 
et al. (1997). 

*^Upper limit on masses based on the detection of 6708 A Li absorption (Kirkpatrick et al. 1999; Reid et al. 
2001 ). 


References. — (1) Close et al. (2002); (2) Lane et al. (2001); (3) Freed, Close, &: Siegler (2002); (4) Reid et 
al. (2001); (5) Koerner et al. (1999); (6) Potter et al. (2002); (7) Goto et al. (2002); (8) Martin, Brandner & 
Basri (1999); (9) Leggett et al. (2001); (10) This paper. 
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2MASS 0559-1404 


2MASS 0937+2931 



2MASS 1047+2124 


2MASS 1217-0311 



2MASS 1225-2739 


2MASS 1237+6526 


Gliese 570D 2MASS 1534-2952 






2MASS 1546-3325 


2MASS 2356-1553 


Fig. Al.— PC chip images around each target source. F814W images are on the left while F1042M images are on the right. Each image 
is 2('3 on a side with a pixel scale of 0(^0455. Image orientations are indicated by the inset arrows, with the arrowhead pointing North and 
orthogonal line pointing East. 







-10 -5 0 5 ±10 -5 0 5 10 

AX (pixels) 


plots of the F814W (top) and F1042M (bottom) PC chip images of 
are 20 pixels on a side, corresponding to 0'/92, and orientations are 
70%, 80%, 90%, and 95% of the source peak are shown. 




F814W 
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0 12 3 4 5 

F814W - F1042M 


Fig. A3.— Color-magnitude diagram of all sources detected at both F814W and F1042M. 3-pixel aperture photometry for point sources and 
target objects (labelled) are indicated by filled symbols, while 5-pixel aperture photometry for extended sources (e.g.., galaxies) are indicated 
by open symbols. The dashed line indicates the limit for detections at both F814W and F1042M. 
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F814W - F1042M 

Fig. A4.— Optical/near-infrared color-color diagram for target objects. 2MASS J-band photometry is used for all target objects. Individual 
photometry for the two components of the 2MASS 1225—2739AB system are derived using Eqn. 2 and the combined light magnitude J = 
15.22di0.05. Data for Gliese 229B are from Golimowski et al. (1998) and Leggett et al. (1999). A straight-line fit to all single sources (excluding 
Gliese 229B) is indicated by the dashed line. 



Fig. A5.— PSF subtraction for 2MASS 1225—2739AB (top) and 2MASS 1534—2952AB (bottom). Both F814W (left) and F1042M (right) 
images are shown. The first image for each set shows the original PG image, while the second shows the residual image after subtracting the 
PSF model. Golor scales are given for each set. 
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Separation (") Seporotlon (") 

Fig. A6.— Completeness limits for companions around 2MASS 0559—1404 in the F814W (left) and F1042M (right) filters. The light grey 
histogram gives the 95% recovery limit; i.e., the limiting flux ratio at which 95% of the simulated binaries were accurately extracted by our 
fitting algorithm at each separation. The 50% recovery limit is shown in dark grey. The separations and flux ratios of 2MASS 1225—2739AB 
and 2MASS 1534—2952AB are indicated by solid circles. 



Fig. A7.— A possible companion to 2MASS 1217—0311. The first two panels show the raw F1042M image frames centered on 2MASS 
1217—0311 prior to cosmic ray correction; the third panel shows the corrected and combined image frame; the fourth panel shows the single 
source PSF-subtracted image, with the putative companion indicated by an arrow; the fifth panel shows the residual image after subtraction 
of two PSFs. All images are 273 on a side, and color scale is indicated on the right. 
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Fig. A8.— Combined mass ratio distribution for T dwarf binaries in this sample and L dwarf binaries in Reid et al. (2001). Individual 
mass ratios are listed in Table 7. The light grey histogram shows the observed distribution normalized to f{q = 1) = 1; number counts are 
shown in parentheses at the bottom of each bin. The observed distribution overlaps a bias-corrected distribution indicated by the dark grey 
histogram and triangles with la uncertainties and upper limits (arrows). The slightly offset filled circles denote the mass ratio distribution 
(with the same normalization) of M dwarf binaries in the 8 pc sample (Reid Sz Gizis 1997), with uncertainties computed as described in the 
Appendix (an upper limit for 0.7 < q < 0.8 is indicated by the downward arrow). The mass ratio limit, qum ~ 0-4, for the combined L and 
T sample for a > 4 — 5 AU is indicated by the short-dashed line. 
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Fig. A9.— Total mass (Mtot) versus separation (a) for star-star (open diamonds), star-brown dwarf (small filled circles) and late-type dwarf 
(primary star later than M8; large filled circles) binaries. The maximum separation for the more massive systems (Mtot ^ 0.4 Mq) appears 
to be limited by log amax oc Mtot (Reid et al. 2001, solid line), while the low-mass binary envelope (dashed line) appears to follow a-max oc 
Mj^j. Resolvable separations for L and T dwarfs (typical distance of 20 and 10 pc, respectively) for 2MASS and HST are indicated along the 
bottom of the figure. 
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Fig. AlO.— Probability distribution for constructed for a sample size N = 10 and number of binaries n = 2. The shaded region gives 
the ±lcr range of acceptable values, whose limits are defined in the Appendix. The integrated probability in this region, 68%, is equivalent 
to l(j Gaussian limits. 



